6 Molecular Orbital Methods 



6.1 Semiempirical methods 

As a general philosophy, the purpose of semiempirical calculations is to 
produce good insight, not to produce bad numbers. 

6.1.1 Extended Hiickel Methods 

Include all valence electrons. Diagonalize an effective, one-electron Hamil- 
tonian 

[H - -EjSjcj = (6.1) 

h ab = (KS ab /2)[h aa + h bb ] (6.2) 
Ka = a a (6.3) 

K is the Wolfsberg-Helmholz constant and has a value of 1.75. 

The overlap matrix S us evaluated using (single- or multiple-^) Slater 
orbitals: 

S ab = (a\b) (6.4) 
The total energy is approximated by the eigenvalue sum 

occ 

E T = Y,E l (6.5) 
i=l 

Carry out Mulliken population analysis. Can be made charge-iterative. 
Historical origins: 

• Mulliken 1949, Wolfsberg and Helmholz 1952; Eq. (6.2) 

• Longuet-Higgins et al. 1954, Lipscomb et al. 1961 

• Hoffmann 1963 

Use in " quasirelativistic mode": 

• Hoffmann school, including band structures. 

For a review on relativistic semiempirical methods, see Pyykko (1988b). 
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6.1.2 Zero Differential Overlap Approximation 

Include all valence electrons. Start from the Hartree-Fock equations 

FVi = [T + V n + V c + V X ]A 

= Etfi + EijiJj (6.6) 

and set systematically all products of two AOs, 

4> A (f>b = 5ab4>l ( 6 - 7 ) 

in the potentials V c and V x , and in S. 
The total energy is calculated as 

E T = (V\F\V) (6.8) 

with the single-determinant wave function 

* = |^i(l)^I(2)... | (6.9) 

Various versions 



• CNDO (Complete Neglect of Differential Overlap): Include only the 
integrals 

(a A (l)b B (l)\r^\cc(2)d D (2)) 

= o~ABO~abO~CDbcd{aa\cc) (6.10) 

• INDO (Intermediate NDO): Add all one-center integrals. Exchange 
introduced. 



Quasirelativistic ZDO models 

• Boca 1987-, R6sch/Zerner et al. 1987-. Use relativistic data. 

• M. J. S. Dewar et al. (1985): Empirical MNDO or AMI parame- 
ters for Hg, Pb. See review by J. J. P. Stewart, J. Computer- Aided 
Mol. Design 4 (1990) 1-105. 

• See Table 7.7 of RTAM III (2000). 

• Thiel and Voityuk (1996): MNDO for d-orbitals, Cl-I, Zn-Hg. 
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6.1.3 Inclusion of Spin-Orbit Splitting 

Many calculations since late 1960'ies, see refs. 69-102 in P. Pyykko, Methods 
in Comp. Chem. 2 (1988) 137-226. Simplest alternative: include the one- 
center SO term, 

centers 

^so = £ C fl l-s (6.11) 

9 

This obviously doubles the Hamiltonian matrix (for a and (3 spin). 
ZDO Methods with SO Splitting 

• Boca (1989, 1990ab): Relativistic CNDO/1. Applications on di- and 
triatomics, Aug + , . . . 

• Bohm et al. (1989): Relativistic ZDO. Fe„, n = 12 - 14. 

• Kotzian et al. (1989, 1991, 1992), Kotzian and Rosch (1991, 1992): 
Relativistic INDO/S - CI ('S' for spectroscopic). Diatomic hydrides 
and oxides. LnO. [Ce(H 2 0) 6 ] 3+ . 

• Jahns et al. (1992): INDO/1 ('1' for ground state) [(Cp^M) 2 C 6 H 4 ], 
M = Sc, Lu. Agostic M. . . H attraction reproduced. 

• Minaev et al. (1989): Relativistic MINDO/3 - CI. Diatomic hydride 
systems. Intensity borrowing. 

• Roszak and Lipiriski (1992): INDO. Triplet - singlet lifetimes of 
pyridine etc. 

6.1.4 Relativistic Extended Hiickel (REX) 

• Introduce \lsjm) basis with separate energy parameters, a a and radial 
parameters d for the SO-split atomic states. 

• Derive these parameters from ab initio atomic calculations, thus ex- 
trapolating them to the molecular domain (Lohr and Pyykko 1979). 
R/NR! 

• Rosch (1983) treated the Kramers degeneracy algebraically (QATREX) 

• Larsson and Pyykko (1986): Charge iterations added (ITEREX) 

• Original default parameters (in 'DEFPAR') acceptable for main-group 
elements. Realistic parameters for the actinides in Pyykko, Laakkonen 
and Tatsumi, Inorg. Chem. (28) (1989) 1801. 

• Always carry out a sensitivity study with respect to parameters or try 
to confirm the ideas found by better methods. 
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Some Results by REX 



• Implement the relativistic theory of nuclear spin-spin coupling con- 
stants. New terms. See Pyykko and Wiesenfield (1981). 

• Symmetry rules for the spin-spin coupling tensor, J, A. D. Bucking- 
ham et al. (1982). 

• Transparent interpretation of the heavy-atom chemical shift in X HX 
NMR, Pyykko, Gorling and R6sch (1987). "HAHA" (1987). 

• The 6p-hole theory of actinyl NQCC, Larsson and Pyykko (1986). 

• Colors of BiPli5 and PbClg~ attributed to the relativistic stabilization 
of the 6s ai LUMO (Schmuck et al. (1990), El-Issa et al. (1991)). 

• The SO splitting of the OSO4 valence 3t2 MO attributed to 5p semicore 
character, in most models. Pyykko, Li, Bastug, Fricke and Kolb 
(1993). 



6.2 One-Electron Molecules 
6.2.1 The Hamiltonian 




(6.12) 



n r 2 



6.2.2 Possible Coordinate Systems 



a) Elliptical: 




(6.13) 



n = -g (S + w 
J2 = ^ (£-v) 



(6.14) 



Constant £: ellipses, constant 77: hyperbolas. 
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Third variable: (p. 



d ^= ?(£ 2 -»7 2 )d£<W (6-14') 



Then 



h^ = E^ (6.15) 
^ = e^/(C,r?) (6.16) 



The Schrodinger equation: 



-m 2 (- r ^- + -^))+^,r ? )-^ /(e,r?) = (6.18) 



1 1 

0. Burrau, Mat /ys. Medd. 7 (1927) 

D. R. Bates, K. Ledsham, A. L. Stewart, Phil. Trans. RSL 246 (1953) 
215. 

See L. Laaksonen et al. (1983a). 
Eq. (6.18) separates, 

f(Z,r,) = L(OM( V ) (6.19) 
The eigenvalues E are obtained from the truncation of the series for 

b) The variable q: 

r)=l-Nlnq, 0<q<l (6.20) 

c) The variable t : 

t = ln£ (6.21) 

d) Prolate spheroidal: 

The 2D-program uses the variables (fJ,,v): 

j £ = cosh/i, < jj, < oo 
[ rj = cos i/, < v < tt 

r\ = ^i?(cosh fi + cos v) 
r2 = -i?(cosh jU — cosz^) 



(6.22) 
(6.23) 
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R 3 

dV = — sinh \x sin ^(cosh 2 fi — cos 2 v)dfidvd(f> 



d 2 



i? 2 (cosh n — cos 2 v) id fi 



d d 2 



+ coth/i- h 



dfi dv 2 



a 2/ i 1 

+ cot^- m — — h 



sinlr/i sini/' 



) SM 



(6.24) 
(6.25) 

(6.26) 



e)Kolb's coordinates (s,t): 

See C. Diisterhoft, L. Yang, D. Heinemann, D. Kolb, Chem. Phys. Lett. 
229 (1994) 667. 



r\ = 4 sinh s + 1 , 



< s < oo 



£ = cost[l + -sin 2 t], < t < 7T 
[Already in L. Yang ei a/., C/iem. P/iys. Lett. 178 (1991) 213] 

6.2.3 Transformation of the Dirac Equation 

D. Sundholm, P. Pyykko, L. Laaksonen, Phys. Scripta 36 (1987) 400. 
D. Sundholm, Chem. Phys. Lett. 223 (1994) 469. 



(6.27) 



Consider the Dirac hamiltonian 



Let 



a ■ p 



V ca-p 
K cq_-pV - 2c 2 

Pz Px ~ ip y 

\Px + iPy -Pz 



The lower equation from (6.28) gives 

ca ■ piph + (V — 2c 2 )V>s = E^ s 
ca ■ p 



2c 2 + E - V 



(6.28) 
(6.29) 



(6.30) 



(6.31) 
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This can be substituted in the upper equation, 

(V - E)ip L = -ca ■ p^s 

= ~ca ■ p(2c 2 + E - V)~ x cq_ ■ p^ h 



Letting 



Recalling 



c 2 



2c 2 + E-V 
h ■ pf(r)a ■ p + V - E]i/j l = 

p = -iV, p 2 = -V 2 , ip = V 



(a ■ A)(a ■ B) = A ■ B + la ■ (A x B) 



QL ■ pf(r)a ■ p = a - [( pf )a- p +fp(a ■ p)] 
A B 

= (pf) -p + ia- (pf x p) + /[ (a-p)(a-p) 

v ' 

p 2 + ia ■ (p x p) 
= 

= fp 2 + (pf) -p + iv- (pf xp) 
= -/V 2 - V/ ■ V - ia ■ (V/ x V) 

In an orthogonal coordinate system, 

wf . IS/,. IS/,. 18/ 
h\ dqi h 2 oq2 h 3 dq% 

(The distance, dsi = hi&qi). Here 

R I 

hfi = K = —y (cosh 2 p, - cos 2 //) 

h<f, = — sinh p, sin v 

where all the hs are weight functions. 
Then 

_ v/ . vv , = _]_ d l d ±_ ]_d£di _ j_d£di 

/i 2 dp, dp, hi dv dv /i 2 dtp dcp 

4 .dfdjj dfdji, 

2 ,. — „2 ..A If),, P.,, Pi,, P), J 



i? 2 (cosh 2 [i - cos 2 v) l d/i dp, dv dv 

4 df dip 

R 2 sinh 2 p, sin 2 v dcj) dcp 
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„ o 4/ ,- coshu d 



i? 2 (cosh 2 /x - cos 2 i/) d^ 2 sinh /j, d/j, 

d 2 cos v d cosh 2 // — cos 2 v d 2 n 
_l 1 1 c 

<9i/ 2 sin 1/ dv sinh 2 ^ sin 2 f d<f> 2 



What, then, is -ia • (V/ x Vip) ? 
On one hand, 



D 



V/ x W = 



df df df 

dx dy dz 

dip dip dip 

dx dy dz 



(all /is equal 1 in cartesian coordinates) 
On the other hand, denoting f x = 



V/ x w 



/t/" ft^^ 
£Vfc ty* 



= e. 



h,,h 



tlebtlu 



+e 4> j^-(f IJi tp u - fvlpp) 
We will need the cartesian components of this determinant, 



D = Di&i + D 2 e 2 + D 3 e 3 
for the scalar product with 



a = o\e\ + (J 2 e2 + 0363 
1 dx dy dz 
1 dx dy dz 

^ = K {ax d^ + ay d^ + az d^ ) 

dy dz 

+ ay d4> + az M ) 



1 dx 

Ha, 06 
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-icr • (V/ X W) = 

dx dy dz 



+ fair. + a yir. + vtjrNMi* ~ fM 



+ fa 



dx 

du 
dx 

8~4> ' " y d(ft 



dy 

' v du 
dy 



Oz—Mf^u ~ fulpfi) 



(6.45) 



From 



x = — sinh /i sin z> cos 6 
2 ^ 

y = — sinh ,u sin i/ sin (6.46) 
R 

z = — cosh fi cos (z is in the molecular axis direction) 
dx R 

— = — cosh a sin u cos 0, ^— = — sinh u cos u cos 0, 

C7/U 2 

dy _ R 

0l ' 2 " " (6.47) 



sinh n sin u sin 
cosh fi sin i/ sin 0, — ^ = — sinh ^ cos f sin 0, 



sinh |U sin cos 
1? 

9/i 2 1 ' dv 2 



dz R . dz R 

sinh \x cos u, — = cosh fi sin u, 



dx 


R 


du 


~ 2 


dx 


R 


d~4> 


~ 2 


dy 


R 


du 


~ 2 


dy 


R 


dcf) 


~ 2 


dz 




du 




dz 


= 




~d~4> 



-^■D = -^-[ (6.48) 

(a x cosh fj, sin u cos (ft + a y cosh \x sin u sin (ft + a z sinh |U cos u) (f u ipcf> — f<j>iftv) 
+{a x sinh fi cos ^ cos (ft + a y sinh \i cos ^ sin (ft — a z cosh ju sin i/) (f^iftn — f^ift^) 
+ {—a x sinh /i sin z/ sin (ft + a y sinh fx sin ^ cos + 0)(f^ipi> — fu^P/j,)] 
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2 ° x 

hfj,h u h^ 



cosh /x sin v cos <\>{Ui>4> ~ U^u) 



+ sinh xx cos u cos <K/</>V7i - //uVty) ~ sinn A* sin ^ sin ~ Mfi) 



„ n — 2 a ?/ 



cosh /x sin z/ sin ~ f$i>v) 



+ sinh cos v sin </>(/^ ~ /juVty) ~ sinh /x sin z/ cos 4>{f^u ~ fvipfi) 



-a z D z = - 



sinh /i cos v{f u ip^ - f^ v ) 



h^h u h, 

- cosh xx sin ^(/^ - f^) 



Recall that 



1 

1 0, 





(6.49) 



-te-D = - 1 -+ T - x 



sinh /x cos v{f v ij}^ - fyipv) 
- cosh /i sin z/^W - /m^v) 



cosh /x sin z/ cos ~ f^v) 

+ sinh /x cos z/ cos (pi f^ip^ - / M Vv>) 
- sinh /x sin u sin 4>{f^ u ~ Mp) 



cosh xx sin z; cos (/»(/^ - /</>Yv) \ 
+ sinh /x cos zv cos 0(/</>Yv - /jiVty) 
- sinh /x sin f sin ~ Mfi) 

-i cosh xx sin f sin 0(/„^, ~ f<$>i>v) 

+ sinh xx cos v sin Yv - f^tp^) 

+ sinh xx sin v cos </>(/;uVv - Mn) 



+i 



- sinh xx cos v{f v ^) ( j ) - /^Vv) 
+ cosh xx sin ~ 7/^) 



cosh /x sin z/ sin <j){f v i>^ - f^u) 
+ sinh n cos i/ sin ^(/^ - 
y + sinh /x sin i/ cos 4>{f^ u ~ Mp) 

It can be noted that the matrix in the above equation is of the form 
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-ia-D = 3 ^ x (6.50) 

Y (cosh 2 fi — cos 2 v) sinh ll sin t> 

/ cosh 11 sin z/(cos ^ - i sin <j)){f v ^ - f^ v )\ 

smh/xcosi/(/^ - ftVv) . . . w . 

, . ,, , , . . +smh ) ucosi/(cos0-ism0)(/^ /1 -/ / xV'0) 
-cosh/zsinz^/^-/^) . , . , . , , . , v , , t , x 

■ smh /i sin z/(sm + l cos (p) {j^ipv ~ JvWn) 



cosh n sin z^(cos ip + i sin 0) (/„*/ty - 
+ sinh n cos z/(cos + i sin 4>) - 



- sinh /i cos v(f v ip<t> ~ f<j>ipu) 




cos0 =F isin0 = ^ , sin0 ± icos0 = < (6.52) 



4i c 2 

-/a • L> = ^ x — — x (6.53) 

R 2 (coslr ji - cos 2 v) sinh ii sin z/ (2<r + E — V ) z 

( a cosh // sin e _1 ^^-V<A ^ 

. . ov du 

smh fj, cos v — ip,!) Qy 

gy — sinh fx cos v -—ip<p 

+ cosh M sin,— ^ r .J S y3 

- sinh /z sin ,(ie --— 

cosh it sin v e l< ^ ^—ips 

ov . , 

- sinh /x cos zv — ip^ 

— cosh /i sin za 

V ^snm M sm nl e ^ — 



— sinh/ucosz/e 1 ^ — — ip^ 

% V9 dVd^ l - cosh/, sin z,^ 

sinh/. sin z, ( ie*)(— --— -U 



. . + cosh// sin - /^) 

\ - smh /x sin z/(sm - l cos (p) - ) / 



Now note that 



c 2 (-i) ay c 2 (-l) 

/m = (2c 2 + J E-F) 2( "^7 ) ' /l/ = (2c 2 + E-V)^~~o^^ (6 ' 51) 

Uoc — = 0, 
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Recalling the Ansatz (6.30), 



/ e i(™-§)<^,L^^\ 

'i(m- ^e^-i^Vi^^) 
i(m+ D^+s^^Ai,!/) 
-(m- l)e^ m -^f(n,u] 



(6.54) 



the operation — \g_ ■ D gives for the upper row: 



a- J r • i 5y i ^ 

Ai< sinh it cos zv— h cosh xx sin v—— 

r , . dV dV 

+ cosh fi sin f — sinh \i cos f — — 

L ov o\x 



e _i *i(m + ^)e i(m+ 5)^2 



sinh /x sin f ie ^e 



V cfyx <9z^ 



A = — 



Ac 2 



i? 2 (cosh 2 /x — cos 2 v) sinh xx sin i^(2c 2 + E — V) 2 



(6.56) 



and for the lower row: 



AU 



. dV dV 
cosh /x sin sinh /x cos i^-^— 



e^ifm 



) e i(m-|)0^L 



+ sinh A ,sin,ie i ^-^(^A-^A)^ 

V Of! OV OV Ofl/ 



dV . dV 

sinh xx cos v— h cosh /x sm z^^r - 



i(m + ^e^+l^ 



(6.57) 
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Dividing away the e^-factors 



-icr • D 



(6.58) 



-A< 



+ 



dV . dV 

sinh a cos v — h cosh u sin z^^— 

. 



(m-^ L 



/ , dV . , 1 

cosh u sin i/— sinh u cos z/—— (m + - ) 

V dv da / 2 



— smh u sin zv — — — 
VO/IO 



z/ dv d\i) J 2 J 



cosh u sin z/— — 
V dv 



sinn /j cos z^— - (m , 



fdV d 
i 



dV d 



+ sinh /j sin zW ^— j- - ^j-) 
\ on dv OV OflS 



dV 



dV 



sinh a cos i/— h cosh /j sin v 

Ov on 



(m + -)V> 



Using the definition of vl, (6.56), 



-io; • £> = + 
/ 



/ 



i? 2 (cosh 2 u - cos 2 v) 2c 2 + E - V 

1, 



"cosf <9F coshu <9F 
+ 



(m 



+ 



shif dv sinh/_i 9// 
vcosh/j dV cosv _|_ 
A sinh u dv sinz^ dfiJ 
fdV d dV d\ 
v d/j, dv dv da) 



/ cosh a dV cos v dV \. 1 . 

I — I ( 777- — — ) 

A sinh n dv sinv dv J 2 
fdV d dV d y 
V dfi dv dv d/i/ - 



\ 



cos v dV cosh a dV 
1 1 

sinz/ dv sinh /j d/j, 



(6.59) 



Combining this result with the diagonal kinetic energy terms (6.38)- 
(6.39) and (V - E)ip, noting that in (6.38) 



-v/ • vv> 



/ 



7 



dV d dV d 



i? 2 (cosh 2 a - cos 2 v) 2c 2 + E - V V d/j, d/j, dv di 
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the upper equation finally becomes 



4/ 



i? 2 (cOsh fjL — cos 2 v) 



+ 



d 2 cosh a d d 2 cos v d 
_l 1 1 1 

sinh/x dfi dv 2 sinz/ dv 



dfi 2 
cosh 2 n 



cos 2 v 



1. 



2c 2 + - y 



sinh 2 /i sin 2 z/ 

/ay a dv d\ 

\ d\x dji dv dv) 



/cosf OF ^ cosh/i <^\/ 
V sin ^ dv sinh fi dfiJ 

/cosh^ ay cost' ay 



2c 2 + £ - y 



V sinh dv 

fdv a ay a \ 



) (m + 
sin ^ a^i / 2 



v a// a^ 
+ (y- 



a^ a/x/ 







and the lower equation becomes 
4/ 



i? 2 (cosh 2 /x — cos 2 v) 



1 



2c 2 + £ - y 



/ cosh /j dV 
V sinh /x a^ 



cos v dV 



(dv a 

V dfi dv 



sin v dfi J 

dv a \ 



) (m — 



2c 2 + £ - y 



az/ a^/ 
/ cosz^ ay coshu ay\ , 

( + -■ -7T + —^ 7r)( 

\ smv ov sinhu oaJ 



m + 



(dv d ay a 

V dfi dfi dv dv 



+ 



a 2 



cosh ad a 2 cos v d 

1 1 1 

a/i 2 sinh/i dfi dv 2 smv dv 

cosh 2 /i — cos 2 v 1 , 

—^2 ~2 ( m + o) 

sinh \i sin v 2 
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